Higher-order photon correlations in pulsed photonic crystal nanolasers 
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We report on the higher-order photon correlations of a high-/? nanolaser under pulsed excitation 
at room temperature. Using a multiplexed four-element superconducting single photon detector we 
measured g (0) with n=2,3,4. All orders of correlation display partially chaotic statistics, even 
at four times the threshold excitation power. We show that this departure from coherence and 
Poisson statistics is due to the quantum fluctuations associated with the small number of dipoles 
and photons involved in the lasing process. 

PACS numbers: 42.55.Tv, 42.50.Lc, 42.50.Ar, 42.55.Ah 



The physics of nanoscale lasers is expected to be rad- 
ically different from that of conventional lasers, due to 
their very high spontaneous emission ratio (/3) into the 
lasing mode, approaching (3 — > 1. Because of the high 
/?, the lasing threshold is low, the system is capable of 
lasing with a small number of emitting dipoles and a 
small number of photons in the cavity mode, and the las- 
ing turn-on time is expected to be very short. However, 
the quantum fluctuations associated with the small num- 
bers of dipoles and photons cause the output intensity 
of the laser to fluctuate strongly with chaotic statistics 
which can be present well above the stimulated emission 
threshold [l|, Q and may even suppress continuous- wave 
lasing Q. Understanding these quantum fluctuations, 
which requires studies of the higher-order photon statis- 
tics of the laser output, is thus crucial for the operation 
of nanolasers, particularly if nanoscale lasers are to be 
used in demanding applications, such as inter-chip opti- 
cal interconnects or integrated clocks. 

Several studies aimed at understanding the physics 
behind the dymics and the noise characteristics of 
nanolasers have been reported in recent years. Measure- 
ments of the second order correlation function g^ 2 '(r), 
which is directly linked to the intensity noise properties, 
were reported some years ago [4J, |5| and more recently 
the evolution of g^ 2 \t, r) was followed with picosecond 
resolution in semiconductor micropillar nanolasers oper- 
ated at 4K in the weak [(| 0] or the strong [1] coupling 
regimes. Lately, the noise properties of nanolasers op- 
erating at room temperature were also reported @. In 
most of these studies, the cavity volume was of the order 
of a cubic wavelength and the gain material consisted 
of a small number of semiconductor quantum dots. In 
some recent work on plasmonic structures, lasers with 
even smaller volumes have been reported [Ig-Qjj], reduc- 
ing further the number of emitting dipoles and photons 
in the cavity. Clearly, in such small lasers, the tradi- 
tional laser theories based on the "thermodynamic limit" 
of large-ensemble averaging and continuum approxima- 



tions for the gain medium and the electromagnetic field 
are no longer valid, as the discrete nature of the num- 
ber of dipoles and photons must be taken into account 
explicitly [Hj]. 

In this Letter we report an experimental study of the 
higher-order photon autocorrelation functions of a single- 
mode photonic crystal nanocavity laser, with quantum 
dots as its gain material, operating at room temperature, 
in the telecommunications wavelength range, and under 
pulsed excitation. In the second-, third- and fourth- 
orders of correlation, the output intensity fluctuations 
display a chaotic component and deviate from standard 
Poisson statistics, even at pump powers of over 4 times 
the threshold. These noise properties can be accounted 
for by the quantum fluctuations arising from the small 
number of photons and dipoles, through a simple model 
that provides new insight into the physics of nanolasers. 
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FIG. 1: Double heterostructure photonic crystal cavity. 
White dashed lines represent the mirror position. The red 
arrow shows the laser cavity 

The nanolaser cavity is formed by a photonic crys- 



tal double heterostructure 
suspended InP membrane 



t] etched on a 320 nm-thick 
and incorporating a single 
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layer of self-assembled InAsP quantum dots [lj| at its 
vertical center plane. The whole structure is grown by 
metalo-organic chemical vapor deposition. The quantum 
dot density is of the order of 1.5 xlO 9 cm -2 and their 
spontaneous emission is centered around 1560 nm at 300 
K with an inhomogeneous linewidth of about 150 nm. 
The cavity is fabricated using electron beam lithography, 
inductively coupled plasma etching and wet etching [16j ■ 
The structure consists of a Wl photonic crystal waveg- 
uide composed of one missing row of holes in the T-K 
direction of a hexagonal lattice structure with lattice con- 
stant a m = 410 nm, and air-hole radius of r = 0.293 a m . 
The lattice constant is modified over two periods at the 
center of the photonic crystal waveguide along the T-K 
direction alone, to the value of a/ = 440 nm, thus form- 
ing a nanocavity with an effective volume of 1.3 (A/n) 3 , 
enclosed by two surrounding "mirror" waveguides with 
smaller lattice constant. The cold cavity quality factor is 
measured to be Q = 45000 [13] , corresponding to a cavity 
lifetime of 1/T C = 35 ps. 

The nanocavity laser is studied at room temperature. 
Optical excitation is provided by a pulsed TkSapphirc 
laser with a center wavelength of 805 nm and a repeti- 
tion rate of 82 MHz. The pump pulse is broadened to 
~ 50 ps after passage through a 25 m length of single- 
mode optical fiber and is then focused to a 5/im spot 
on the sample with a microscope objective (NA = 0.4). 
When the quantum dots are highly excited and contain 
several electron- hole pairs, they emit into a broad spec- 
trum with a lifetime of I/711 = 400 ps 0, [l8| which 
feeds the cavity mode (l9^ . Emission from the nanocav- 
ity laser is collected with the same objective, directed 
through a dichroic beamsplitter that blocks the pump 
light and sent to a 0.75 m grating monochromator. One 
output port of the monochromator holds an image in- 
tensified near-infrared camera, which records the emit- 
ted spectrum and power. Alternatively, an internal mir- 
ror can be flipped, directing light to the other output 
port, where it is filtered with a slit and subsequently 
collected into a single-mode optical fiber and sent to 
the detector. The detector is a four-element supercon- 
ducting nanowire single-photon detector (SNSPD) j^j 
in which four independent, single-photon-sensitive ele- 
ments are interleaved over a single spatial mode of the 
optical beam [21] . Each element consists of a current- 
biased superconducting nanowire that is driven into a 
resistive state, thus delivering an output voltage pulse 
upon absorption of one photon. The active area of the 
four-element SNSPD is matched to the mode diameter 
of the fiber, allowing all four interleaved nanowire ele- 
ments to equally sample a single spatial mode [2jJ, [22 1 . 



Fast four-channel electronics record photon arrival times 
on each element. These time-tag data are post-processed 
to obtain multi-start, multi-stop correlation histograms 
between two, three and four SNSPD elements for a con- 
tinuous range of all associated time delays. The SNSPD 



is held at a temperature of «3 K in a closed-cycle he- 
lium cryocooler. Each element has a system detection 
efficiency of » 1% at 1550 nm and a dark count rate of 
~ 100 Hz. Details on the data processing, including the 
elimination of crosstalk among the interleaved detectors, 
are given in the annex on supplementary information [23| . 
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FIG. 2: Third order autocorrelation function g' 3 ' (n, T2) at 
a pump power P=2 Pt^- The tallest bar is located at ri = 
t-2 — and represents <?' 3 '(0, 0). The three ridges of elevated 
bars that fall along n = 0, T2 = and n = T2 have values 
approximately equal to (?' 2 '(0). The remainder of the bars 
have an average value of ~1. 

Figure [2] presents a typical set of normalized exper- 
imental three-photon coincidence data, extracted from 
the time-tagged arrivals of photons on the four multi- 
plexed SSPD detectors. It corresponds to the third order 
autocorrelation function, g^ 3 \ri, T2), in histogram form. 
The blocks at ^ t\ 7^ t 2 7^ are normalized to an aver- 
age value of one as the three photons detected belong to 
different excitation cycles; the blocks at T\ = 0, at r 2 = 0, 
or at T\ = t 2 correspond to two photons being detected 
in the same excitation cycle (and the third one in an- 
other cycle) and thus give directly the value of </ 2 )(0); 
finally, the block at t\ — t 2 — 0, where all three photons 
are detected in the same excitation cycle, gives the value 
of 3 (3) (0,0) [53. For fully coherent light g {n) (6) = 1, 
whereas for chaotic (thermal) light g^(0) = nl, for all 
orders n. For partially coherent light, we may thus define 



g {n) (Q)-i 
n\-l 



(1) 



which is a normalized ratio, indicating if the statistics 
are near coherent or chaotic light : Below threshold 
the output light is chaotic, so that h^ n ' — 1, whereas 
above threshold, one expects a coherent laser output and 
-> for all n. 
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On Figure [3] we present the experimental data and 
the theoretical curves (obtained below) for the Light- 
in Light-out curve (right hand axis) and ftW for n = 
{2, 3, 4} (left hand axis) as a function of the pump power 
of the nanoscale laser under pulsed excitation. The lasing 
threshold is clearly identified at the "knee" of the Light- 
In Light-Out curve (which corresponds to the inflection 
point of the log-log curve), at P t /j=100 [iW, while the (3 
factor is given by the ratio of the two slopes of the Light- 
In Light-Out curve (corresponding to the "jump" of that 
curve m log-log scale) of /3 = 0.008. The ftW, which 
below threshold are expected to be equal to 1, steadily 
decrease to reach = 0.05 ± 0.04, = 0.036 ± 0.02 
and = 0.029 ± 0.015 at P = AP th = 400^W, all 
significantly higher than the value of expected for co- 
herent emission. The corresponding g^ n '(0) are g^ 2 \0) = 
1.05±0.04, 5 ( 3 )(0) = 1.18±0.10,g< 4 )(0) = 1.67±0.34. In 
other words, the intensity fluctuations still have a chaotic 
component at four times the threshold power, in sharp 
contrast with the behavior of conventional lasers (with 
/3 « 10~ 5 ), for which all rapidly vanish above thresh- 
old. 



number of photons in the cavity, s, be integers. 
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FIG. 3: (Right axis) Experimental (full dots) and calculated 
Light-In Light-Out curve. (Left axis) Experimental values of 
h (n) (0) as function of pump power. Open circles, squares and 
triangles are for n=2,3,4 respectively. Full lines represent the 
expected ^"'(O) without any free parameters. 



The persistence of the chaotic component can be sim- 
ply understood following Refs. Because of the rela- 
tively small numbers of the emitting dipoles and photons 
in the cavity, a small fluctuation in these numbers in- 
troduces a non-negligible change in the stimulated emis- 
sion rate, inducing chaotic (thermal) statistics even above 
threshold. 

More quantitatively, we may model our system in 
terms of the traditional laser rate equations [2J under the 
constraint that the number of excited dipoles, N, and the 



dN(t) 
dt 

ds(t) 
dt 



= Pf(t)- 7{l N(t) 
= -T c s(t) + /3 ll{ N(t) 



-p 7ll N(t)s(t) (2) 

-frnN(t)8(t) (3) 



where f(t) is the temporal profile of the pump pulse in- 
tensity, of integral equal to one. For each incident power 
P, the rate equations are solved numerically by iterat- 
ing over a time-step dt which is small compared with the 
characteristic times of the system and l/r c , and us- 
ing a quantum jump approach 24| to ensure that TV and 
s are integers: At each iteration step, the right hand side 
of Equations ([2} and (|3]) is evaluated, taking into account 
the discrete evolution of N and s through a sequence of 
three binomial processes. First, the number of dipoles 
having decayed between t and t + dt is drawn from a bi- 
nomial distribution with probability p\ = 711 (1 + f3s(t))dt 
for each one of them to decay. Second, the decaying 
dipoles each produce a photon that has a probability 
P2 = ^i+ s ps(t) °f entering the cavity. And third, the 
photons that have accumulated in the cavity can escape 
with probability = T c dt. Note that for consistency, we 
also take into account Poissonian statistics for the pump, 
although it does not alter the overall conclusions. In the 
end, for each incident power P, the g (n) (0) are calculated, 
according to 



ff (n) (0) 



</ t nto" 



\ Si {t)-k) > 



(</**(*) >)» 



(4) 



where Sj (t) is the i — th realization of the number of pho- 
tons in the cavity, arising from the rate equations, and 
< x > denotes the average of x over all realizations. Av- 
erages over typically 1000 realizations are performed. As 
can be seen on Fig. |3l the experimental data are very 
well described by our model for the Light-In Light-Out 
curve as well the second and third order correlation func- 
tions. They are reasonably well described for the fourth 
order correlation function, considering the very large es- 
timated experimental error associated with the measure- 
ments of g( 4 ) (0). Details of the calculations can be found 
elsewere [25|. Here, it is sufficient to point out the main 
differences with other methods for calculating the (0) . 
For traditional lasers with a large number of dipoles and 
photons in the cavity, quantum noise is introduced in 
the classical rate equations by means of a Langevin driv- 
ing force with Gaussian statistics 0, Q] . This approach, 
however, does not take into account the discrete nature 
of the particles and is thus not suited for lasers with a 
small number of dipoles and photons. In addition, the 
usual Langevin force statistics have only a second cumu- 
lant. This means that the higher-order averages can only 
be evaluated under specific assumptions, for example, if 
they correspond to sums of products of lower-order aver- 
261 ] . For small lasers, the discretization of the num- 



ages 

bers of dipoles and photons is usually treated through a 
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set of master equations for the probabilities of the states 
with a g iven number of excited dipoles N and photons 
s 3, 13]. While this approach is ideally suited for very 
small systems, it requires calculating N§ interdependent 
probabilities, while our method requires the computation 
of trajectories (realizations) with at most Nq jumps. Sec- 
ondly this approach allows for a more intuitive approach 
to understanding the underlying physical phenomena. 




300 400 500 600 700 800 
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FIG. 4: Dashed lines: 4 different realizations of the stochastic 
rate equations at P — l.SPth- Full line: The mean pulse 
averaged over 1000 realizations. It can be clearly seen that 
the longer it takes for the pulse to emerge, the smaller its 
amplitude. (Inset) Dots: Statistical distribution of the total 
number of photons per pulse at P — 1.8Pth- Line: Expected 
distribution for Poisson statistics. 



wave lasing due to fluctuations in small lasers, discussed 
by Roy-Choudhoury et al. [3|. Fluctuations may sup- 
press a nascent laser pulse, leading to "false starts" and 
delaying the time when the successful laser pulse takes 
off. This asymmetry in the statistics of pulse intensity 
can be characterized by the higher (third and fourth) 
moments of the distribution and thus impacts the higher- 
order autocorrelation functions. 

In conclusion we measured up the the fourth order au- 
tocorrelation function for a photonic crystal nanolaser 
operating at room temperature under pulsed excitation. 
We observe that all values of the autocorrelation function 
are significantly above unity even at 4 times threshold, 
indicating the presence of chaotic fluctuations in spite 
of the predominance of stimulated emission. A simple 
model taking into account the discrete nature of the num- 
ber of photons and dipoles succesfully describes the ex- 
perimental data, providing a novel understanding of the 
operation of the nanolasers. Measurement of the higher- 
order photon correlation functions thus opens the way to 
the study of very small lasers, whose operation deviates 
strongly from the "thermodynamic limit" that governs 
conventional lasers. 

The authors acknowledge financial support from the 
Triangle de la Physique under the BIRD project and from 
the French National Research Agency (ANR) through 
the Nanoscience and Nanotechnology Program (project 
NATIF ANR-09-NANO-P103-36). 



To gain a better insight into the effects of fluctuations, 
we present in Figure HI the evolution of the number of 
photons s(t) for four realizations of the rate equations 
(dashed lines) at an excitation power P = l.8Pth as 
well as the mean of s(t) over 1000 realizations (full line). 
We observe that even at pump power two times that of 
threshold, there are strong intensity and timing (jitter) 
fluctuations, which are correlated: The longer it takes 
for the pulse to be formed, the smaller its intensity, be- 
cause in the meantime the number of excited dipoles has 
decreased. The inset of Figure 0] represents the statisti- 
cal distribution of the integrated pulse area (in number 
of photons), fitted by a Poisson distribution of the same 
mean value. Clearly the number of photons in a pulse 
deviates from a Poisson distribution, and this will cause 
the autocorrelation functions g'"' (0) to deviate from the 
value of one characteristic of coherent emission. As can 
be seen in the inset, the deviation from a Poisson distri- 
bution consists of a lower probability for high-intensity 
events and a higher probability for low-intensity events. 
The correlation between the intensity fluctuations and 
the jitter implies that late starts are more probable than 
early starts, even though the gain medium can provide 
more photons at earlier times. The physical origin of this 
effect may be similar to the suppression of continuous- 
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